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Abstract 

The problem of oscillating flows inside pipes under periodic forcing of viscoelastic fluids is ad- 
h— 5 dressed here. Starting from the linear Oldroyd-B model, a generalized Darcy's law is obtained in 

qq frequency domain and an explicit expression for the dependence of the dynamic permeability on fluid 

parameters and forcing frequency is derived. Previous results in both viscoelastic and Newtonian 
fluids are here shown to be particular cases of our results. On the basis of our calculations, a possible 
explanation for the observed damping of local dynamic response as the forcing frequency increases is 
given. Good fitting with recent experimental studies of wave propagation in viscoelastic media is here 
exhibited. Sound wave propagation in viscoelastic media flowing inside straight pipes is investigated. 
In particular, we obtain the local dynamic response for weakly compressible flows. 

^ 1 Introduction 

o 

Wall-bounded oscillating flows of Newtonian fluids driven by periodic forcing have been extensively stud- 
ied both theoretically and experimentally. Interesting phenomena such us the formation of non-steady 
^Ch boundary layers have been observed, provided that the frequency of the applied force is larger than the 

inverse of the characteristic time in which the vorticity diffuses along the cross-section (see for instance [1] 
and references therein). Recently, the formation of a non-steady boundary layer has also been predicted 
t-H for Oldroyd-B fluids g]. 

In recent advances in many topics ranging from soft-matter, biofluid mechanics, rocket propulsion, 
rheology among others are encountered flows in which Newtonian approximation breaks down. Remark- 
ably, it has been recently observed in oscillating wall-bounded flows of these fluids that when relaxation 
times of elastic stresses are comparable with characteristics times scales of the flow, the nontrivial inter- 
actions of elastic degrees of freedom with the inertial ones and flow geometry give rise to a resonant like 
behavior that is absent in the corresponding Newtonian flow ([3],[1],[5]). Thus, the description of oscil- 
(f-) lating wall-bounded flows of non-Newtonian fluids is of great importance for both fundamental physics 

t— I and applied topics. 

In Refs. (|1],[5]), by means of velocity measurements at the center of a viscoelastic fluid column 
(NayCl/NaSal in water) flowing inside a straight pipe of circular cross section [3] and more recently 
in the entire cross section of the fluid column [5], it was found a dramatic enhancement in the dynamic 
response to an oscillating periodic pressure gradient, that is, the fluid response, measured in terms of 
the velocity for a given amplitude of the pressure gradient, exceeds by several orders in magnitude that 
obtained for the steady case. It is noteworthy that most of the features observed in experiments, namely, 
the values of the driving frequencies at which the resonant behavior takes place and the formation of 
Couette-like flows as the driving frequency increases are successfully reproduced by a linear Navier-Stokes 
and Maxwell constitutive equations ([3], [5], [5]). Nevertheless, theoretical predictions within this model 
overestimate the magnitude of the response amplitude as well as the root mean square velocity at the 
pipe axis and the velocity at certain points of the cross section of the fluid column. Furthermore, the 
Maxwell model cannot describe the damping observed in the peaks of dynamic permeability as the driving 
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frequency is increased (at least in the range of frequencies covered by both experiments). These issues 
motivated us to extend earlier theoretical works in the hope that this may shed light on the the physical 
mechanisms that could lead to the observed phenomena. 

In this article, specifically, we shall try to assess the effects that pure Newtonian contributions (due 
to the solvent) to the total stress tensor of viscoelastic solutions will have on the fluid response in wall- 
bounded oscillating flows under periodic forcing. In order to achieve that, we shall assume that stresses 
in the viscoelastic solution obey the linear Oldroyd-B constitutive equation ([7], [E],[S])- In particular, 
we present a detailed theoretical derivation of the velocity field and the dynamic permeability in flows at 
small both Reynolds and Weissenberg numbers. As we shall see, our finding suggests that the (additional) 
Newtonian component of stresses due to solvent viscosity might leads to the damping of the local dynamic 
response as the forcing frequency increases, a fact observed in ([4], [5]). We also investigate the behavior of 
these quantities in weakly compressible flows. Our main purpose is compare our findings against known 
results for incompressible fluids. 

The article is organized as follows. In section II, for completeness, we briefly review main aspects of 
the Oldroyd-B model. In section III we introduce the main assumptions and the resulting set of governing 
equations we have used to derive explicit expressions for the dynamic permeability, the velocity field and 
the local dynamic response which are defined within this section. Here we also compare our results with 
the measured values of the local dynamic response and the root mean square velocity reported in 0]. In 
section IV, the problem of modeling stresses for weakly compressible flows is revisited and starting from 
it, a explicit expression for the local dynamic response is obtained. A comparison of the latter magnitude 
against the corresponding incompressible case, bearing in mind the experimental conditions, is also made. 

Some of the definitions and expressions presented in this article have already been presented in previous 
publications. Wherever needed, we shall repeat them for the sake of clarity and completeness. 

2 Oldroyd-B model 

One of the simplest models that describes basic viscoelastic behavior of an incompressible liquid having a 
constant shear viscosity, and constant relaxation and retardation times is the so-called Oldroyd-B model 
(0, [H] 5 ) - This is an extension of the Maxwell model that fits quite well data from polymeric solutions. 

In the Oldroyd-B model, it is considered that viscoelastic fluids can be regarded as a dilute suspension 
of elastic dumbbells. Thus, the complicated structure of the molecules that form the liquid is replaced 
by a simpler mechanical structure consisting of identical pairs of microscopic beads (of negligible mass) 
connected by Hookean springs. Besides, the concentration is supposed to be uniform and very low. The 
latter ensures that polymer-polymer interactions are negligible. Furthermore, the viscoelastic solution is 
regarded as a continuous medium. The "coarsening" involved in such description requires to get rid of 
the microscopic degrees of freedom. Consequently, the reaction of the dumbbells on the fluid is treated 
at a mean field level described by an elastic contribution a e to the total stress tensor, which is found to 
be proportional to the conformation tensor, i.e. 

n\, x RiRj >, (1) 

where Ri denotes the i-th component of the elongation vector (obtained by subtracting the position of 
one bead from that of the other one). In equation [I] the average is taken over the statistics of the thermal 
noise, or equivalently over a volume large enough to contain a huge number of molecules but small 
compared with the whole fluid volume. The conformation tensor is defined as of - = — ' 2 J , where R is 
the equilibrium length of the spring. The evolution equation for the conformation tensor can be inferred 
by considering the forces which act on the beads, namely, the hydrodynamic drag which obeys a Stokes' 
law, a Brownian force and the elastic spring force. When details of the kinetic theory are worked out one 
gets the following model incorporating the elastic nature of the dumbbells stress tensor ( [5] , [TU] , [TT] ) , the 
so-called Maxwell model 



a e + = Vp [Vw + (V«) 4 ] , 



(2) 
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where rj p is the contribution of the elastic additives to the total shear viscosity in steady flows at small 
shear rates, (Vw)ij = f^r is the velocity gradient tensor, (Vw) 4 is its transpose and A is the relaxation 

time of elastic stresses. The symbol 4r denotes the upper-convected time derivative defined by ([9"].|12|) 

lt =^ t +(W)G-[G-Vv + (VvY ■ G] (3) 

Thus, the above time derivative takes into account that the relation between stresses and kinematic 
tensors at a fluid particle are independent of the instantaneous orientation of that particle. One can 
see that the ratio between the terms inside brackets in the r.h.s of equation [3] to the linear relaxation 
term is determined by the dimensionless expression Ay, called the Weissenberg number, where V and L 
are a typical velocity and length scale of the flow, respectively. For instance, in a flow inside a straight 
pipe those magnitudes could be taken as the maximum value of the velocity across the cross section and 
the radius, respectively. We may also note that the ratio between the advection term (v ■ V)G and the 
relaxation one is given by the Weissenberg number as well. Thus, when the latter number is sufficiently 
small the relaxation of elastic stresses overcomes the effects caused by nonlinearities in the constitutive 
equation [Eqs. [2] and [3j. In this limit, the upper-convected derivative could be replaced by the local time 
derivative 

- - -• (4) 
St dt { ' 

In the Oldroyd-B model the total stress tensor a is given by the sum of a Newtonian solvent contri- 
bution er s and the elastic additives contribution a e , i.e. 

a = (j e +(T s . (5) 
The constitutive equation for the solvent is given by 

a s = Vs [Vw + (V«)*] , (6) 

where rj s is the solvent viscosity. 

When the two contributions are added, the result is found to be 

( D+ *££)- m 

where D = Vw + (Vi>) is the symmetric part of the velocity gradient tensor, and t)q = rj s + rj p . 

From the above expression the Maxwell's model is recovered making rj s — > 0. Moreover, if we let 
A —> 0, the Newtonian model is obtained. 



3 Non steady flow of Oldroyd-B liquids in straight cylindrical 
pipes. 

3.1 Set of governing equations 

We shall now proceed to derive an analytical expression of the velocity field of an Oldroyd-B liquid subject 
to an oscillatory pressure gradient and confined in a pipe of uniform circular cross section. 

In order to achieve it we must solve the Navier-Stokes equation for incompressible fluids, i.e. 

p(^ + (v-V)ti\ =-Vp + V-a (8) 

V • v = (9) 



with a given by the constitutive equation 12 This is a highly nonlinear set of coupled equations. Hence, 
it turns out to be rather difficult to carry out an analytical treatment of any flow property of a viscoelastic 
fluid. 
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Henceforth, we shall consider only flows in which the Weissenberg and Reynolds numbers are very 
small. The latter means that the convective term (v ■ V)w in the Navicr-Stokcs equation is negligible. 
With these assumptions the set of governing equations is determined by 

dv 

p— = -Vp + V-a, (10) 

V-u = (11) 

and 

Hit)- < i2 > 

A more precise discussion of the above assumptions comes in the following. This type of flow occurs, 
typically, under the influence of a reciprocating piston located at one end of a pipe of circular cross- 
section ([IJ,[4], [5] and references therein). We assume that pipe length is much larger than its radius 
and we shall study the flow very far from the moving piston. There, under the previous assumptions, 
the flow can be described by the component of the velocity along the pipe axis. Moreover, owing to the 
incompressibility assumed and symmetries involved, the velocity only depends on the distance to the pipe 
axis, v = v z (r,t)e z and the pressure only varies in the longitudinal direction Vp — dp ^'^ e z . With these 
assumptions the pressure gradient must be constant along the pipe. 

Then, the convective term in the Navier-Stokes equation vanishes identically leading to 

,*^ = -Vp + V-a. (13) 



Taking divergence in both sides of equation (12), from the definition of D, and bearing in mind the 
expression for the velocity, it can be obtained 2V • D = V 2 w z . These assumptions lead equation 12 



to 



V^ + A— = ),«, V\ + A r -^|, (14) 

where A r = A— is a constant with dimensions of time. 

If we derive with respect to time, multiply by A both sides of ( 13 1 and substituting the term A 



from ( 14 ) we obtain 



dt 



d 2 v z (r,t) dv z (r,t) OVp ( 2 . 9V\, 

X P ^^+P = - X ^ +r lo V v z + \ r —t—\ (15) 



OP 1 ^ dt dt ' lu \ z ' ' dt 

Let us now refer all linear dimensions to the pipe radius R, time to — and the velocity to where K 
is the amplitude of the pressure gradient divided by the fluid density and v Q is its kinematic viscosity at 
steady motion. 

Thus, it can be obtained 

„<9 2 A(£,r) <9A(£,t) - _ N A VsdV 2 A ( A dVp _ A 

where A = ^j 1 denotes the inverse of the Deborah's number, r, £ and A(£, r) denote the dimensionless 
time, radial coordinate and velocity, respectively. 

Fourier transformation of the latter equation leads to 

1 d 0A(f,a) 1-iaA ~ 1-iaA dp(a) 

kz — ) + ia- — — A(^,a) = - —r-,—^. — (17) 

4 at, at, 1 — iaAr]' 1 — xaArf oz 

where a = This number is the ratio between inertia and viscosity forces, so it can be regarded as 

a Reynolds number for non steady flows. A(£, a) and denote the Fourier transform of the velocity 

and pressure gradient, respectively. Hereafter, we shall drop the upper hat and the dependence on a will 
indicate the Fourier Transform, unless otherwise noted. In equation ( j 1 7j ) r\ — 2& is the viscosity ratio. 



To obtain the velocity field, boundary conditions must be 
the fact that the velocity must remain finite at the pipe axis. 
Then it follows for the velocity 
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added, namely, the no slip at the wall and 



H ^) = -[^-^ m ] j^r> ( 18 ) 

't , , 1 — iaA 
/3 2 (a)=ia- — - 19 

If we make r\ — > in the above equation the same result as for the Maxwell model is obtained, (see 
Refs. [[3] ,(4] , [E]]). In addition, making also A — > the Newtonian behavior is recovered, (see p]). 

Now, following the same procedure as in Refs. ([3], [3], [5]) we shall calculate the average velocity over 
a cross section. Thus, 

where the symbol ((■■■))£ denotes the spatial average over the cross-section of the quantity inside brackets. 

From the above expression it can be inferred that the total flux is proportional to the pressure gradient. 
This resembles Darcy's Law in frequency domain. Then, we define the dynamic permeability as 

3.2 Comparison with experimental results 

To compare with experimental results we shall, as in [4], define a local dynamic response (LDR) as 
follows 

/ x P 2 A(0,q) 
K(a) = -R 22) 
dp{a)/az 

In Ref. [3] a detailed experimental study of the dynamic response of a Newtonian fluid, the Glycerol, 
and a viscoelastic fluid, CPyCl/NaSal solution, under an oscillating pressure gradient was performed. 
Measurements of fluid particles velocity and the root mean square velocity V rms were made inside a 
straight vertical cylinder with circular cross section. 

In Ref. [1] the pressure gradient was described by 

dp(t) 2 



dz 



pzow sin{u)t) (23) 



in which zq represents the piston displacement amplitude equal to 0.8mm. Thus the experimental value 
of the (LDR) K e (a) is defined in Ref. [3] as 

K ( a ) =r lo-, tt- (24) 

dp rms /dz 

Hence, since the velocity in the experiment varies sinusoidally, the above expression assumes the form 

(25) 

where v denotes the amplitude of the velocity oscillation. 



Thus, the absolute value of (22) can be compared directly with (25), provided the velocity and the 
pressure gradient are sinusoidals. 

Theoretical results given by the Maxwell model ( [4] , [5] , [13] ) , reproduce quite well the values of a at 
which the peaks in the LDR of the CPyCl/NaSal solution are observed. However, the predicted relative 
amplitude deviates from the experimental results. In the Maxwell model the relative amplitude is larger 
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than the experimental values and this model also fails in predicting the decrease of the LDR observed 
experimentally. In Ref. ?j it is speculated that these differences might be due to nonlinear phenomena 
that occur in the CPyCl/NaSal solution as well as to compressibility effects that could take place when 
the frequency is increased. Regarding the nonlincaritics, as it was mentioned in the Introduction, there 
are two potential sources of them in the hydrodynamic equations, inertial and elastic. The former source 
is vanishingly small, since the Reynolds number is set R e <C 10 -4 (see [3]). As for the elastic sources 
of nonlinearity we must also say that, in [5], measurements were carried out at a more detailed level, 
namely, the whole velocity field was explored along the cross-section and at different heights from the 
moving piston. The former measurements confirm the hypothesis (to a certain extent) that the velocity 
field depends only on the radial coordinate. With this assertion in mind, let us now examine a little more 
closely the Oldroyd-B constitutive equation (Eq. [7]) in the case of pure shear flows. 

Let us first note that in such flows the velocity gradient tensor satisfies, besides the incompressibility 
condition TrD = 0, an additional constraint (V») i3 - (Vw) jfe = [2]. Moreover, the operator v-V vanishes 
identically, as mentioned above. Then, it follows for the upper convected derivative of the velocity gradient 
tensor 



5D _ dD 

With the latter result, the Oldroyd-B constitutive equation can be written as 



2(V«)' ■ Vv. 



r)rr 

X— - A [<j • V« + (Vv)* • a] 



D + A 



dD 
~dt 



2(Vv)* • Vu 



(26) 



(27) 



The effects (assumed weaks) of the nonlinear terms can be assessed by a perturbation procedure. 
Indeed, an expansion of the stress tensor in terms of the Weissenberg numbered = <To + (Wi)ai+0(Wi 2 ) 
substituted in Eq. [27| gives to zero and first order 



o-q + A 



do-p 

dt 



2t7o(Vw + (V«)') 



and 



o-i + A 



~dt 



0( W + (V*V)«) 
dt 



2»j.7 V*«* ■ (W*v*Y + [(T ■ W + (W)' ■ <7 ] 



(28) 



(29) 



respectively. In equation |29[ the quantity V*u* denotes a dimensionless velocity gradient tensor defined 
by the relation Vw = 7 V*i>*. In the particular case we are concer ned wit h, in which the only non- 



f^from Eqs. 



28 



and 



29 



we can see that the last 
stress a zz . This shear-induced 



vanishing component of the velocity gradient is Vi>, 

two terms in the r.h.s of [29] are responsible for the appearance of a norma 
anisotropy in the stress tensor is a nonlinear effect which is absent from the corresponding Newtonian 
flow and has been disregarded in our present derivations. Due to this anisotropy, if the fluid suffered 
a perturbation of large enough amplitude it would have a nonlinear (subcritical) transition to a sort of 
weakly turbulent state |15j . In such a case our analytical treatment will break down. However, from the 
analysis of the order of magnitudes of the terms involved in the equations resulting from our perturbation 
procedure, we can note that shear stresses becomes much larger than the normal ones when the following 
condition is satisfied ^ < 1. Then, the replacement of the upper-convected time derivative by the local 
time derivative seems to be a good approximation at high frequencies and indeed, as we shall see, our 
linear theory provides a better fit at high frequencies. Nevertheless, at low frequencies non linear effects 
must be taken into account and in our opinion, are the main cause of the evident disagreement with the 
predictions of both the linear Oldroyd-B and Maxwell models, as we shall show in the following. 

In Figure 1 the comparison of our calculation with experimental results is shown. As in Ref. [4] 
we plot the LDR scaled by its value at steady motion versus a. The theoretical curve exhibits a good 
agreement with experiments. The curve shape is quite similar for values of a in the range of 0.4 to 1 and 
only slightly differs in the amplitude of the response. We can see that the decrease in the dimensionless 
LDR observed experimentally is quite well reproduced. This allows us to infer that as a increases, the 



1 \n oscillating flows we can define the Weissenberg number as Wi = \"io, where -f is the root mean square of the 
characteristic velocity gradient 
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purely Newtonian component becomes stronger and dissipative effects are more evident. In fact, this can 
be predicted from equation (12 1. From the definition of a and assuming that the viscosity r\ does not 
depend on the velocity gradients, an increase in a leads to an increase in the number wA^J. If the latter 



number is much larger than unity we have in equation ( 12 ) 

da 



» * 



IIA 



dD, 



» IIX>I 



dt rjo dt 

Therefore, in this limit we have for the stress in the viscoelastic solution 

a ~ 2r) s D 



(30) 



(31) 



Hence, the solution behaves in this Newtonian liquid. Moreover, from the above arguments 

we see that the Newtonian factor becomes stronger as the number wA— is increased. This is a significant 
difference from Maxwell liquids, whose behavior is similar to a Hookean solid when a; A 3> 1 [9j. Thus, there 
is a transition from viscoelastic to purely dissipative (Newtonian) dynamics of oscillating flows of Oldroyd- 
B fluids. We stress the fact that the relevant parameter for this transition is the number wA 2 ^, which does 
not depend on the geometrical properties in which the fluid is confined. Another noteworthy observation 
is that, neither nonlinearities nor compressibility effects have been taken into account. Nevertheless, 
in Figure 1 it can be seen that both models exhibit similar behavior and overestimate the relative 
amplitude of the LDR at low frequencies. These disagreements, as mentioned above, might be caused by 
nonlinearities. 

The curve for the Oldroyd-B model was obtained using the relaxation time, A = 1.9s, and the viscosity 
at steady flow, 770 = 6QPa ■ s reported in Refs. ([I], [IE])- Besides, the value of the solvent viscosity which 
provides the best fit to the experimental data was i] s = 0.08Pa ■ s. This value is approximately 80 
times larger than the water viscosity, which is the solvent of the CPyCl/NaSal 60/100 solution. Thus, 
the Newtonian factor considered here within the Oldroyd-B model must be thought as a contribution 
leading to an effective viscosity larger than the water viscosity. This result might be due to the fact 
that CPyCl/NaSal 60/100 solution is not very dilute jTHj, whereas the Oldroyd-B model is valid only for 
dilute polymer solutions ([12], [17]), where the hydrodynamic interactions are neglected. Nevertheless, 
for highly concentrated solutions, disturbances of the solvent velocity field generated by the motion of the 
polymer molecule (coil) affects the drag force on the neighboring coils. As a consequence, the total drag 
force is larger than the total force that would be produced if the solution were very dilute. So, using in 
this case the Oldroyd-B model leads to consider effective values of the parameters involved in the model. 

1000 
500 
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Figure 1: Comparison between the Maxwell and Oldroyd-B models with experiments for the CPyCl/NaSal 
solution . The dashed line is the absolute value of k(o) predicted by the Maxwell model. The continuous 
black line is the result for the Oldroyd-B model with the same value of A=173.71 and ry = 60Pa ■ s. The 
solvent viscosity used was ry s = 0.08Pa • s, leading to i] = 0.8 • 10~ 3 . Experimental values are shown by 
points. 

Another quantity reported in Ref. [4] is the root mean square value of the velocity at the pipe axis.. 
Since, as mentioned before, the pressure gradient varies sinusoidally we assume that the velocity field has 



the form 

v z (r,t) = Im((f>(r)e iut ) (32) 

Substitution of equations (32) and (231 in (15 1 leads to an ordinary differential equation for </>(rFl 
whose vanishing solution at the pipe wall, r = R, is given by 



MPr) 



ZqUJ 



1 = -i ' 



1 + iuj\ 
v V 1 + iu\r/ s /r)o 



From equations ( 33 ) and ( 32 ) we have for the velocity 

v z {r,t) = Ho(r,u])sin(£ + cot) 

where H (r : oj) = ^ Re 2 {4>{r)) + Im 2 ((j)(r)) and ( — ArcTan j^^fey] 
Therefore, the root mean square of the velocity can be calculated as 



(33) 
(34) 

(35) 



V„ 



V2 



(36) 



The next figure shows a comparison between predictions of the Oldroyd-B model and experimental 
values for the root mean square velocity at r = 



~: o.io 




Figure 2: Superposition of the root mean square velocity predicted by the Oldroyd-B model at the pipe 



axis represented by the continuous black line, equation (36 1 evaluated at r = 0, with the experimental 
values, shown by points in the graph. A good fit is here exhibited for frequencies higher than 4 Hz. Peaks 
are quite well predicted and only differences are observed in their magnitude. The values of density and 
viscosity of the CPyCl/NaSal solution were used ([3], [IS]). The solvent viscosity that provides the fit is 
equal to 0.08 Pa ■ s 

Figure 2 shows that the Oldroyd-B model can reproduce accurately the frequencies at which peaks of 
the root mean square velocity are observed. Moreover, the curve shape is well reproduced for frequencies 
above 4 Hz. These agreements reinforce the fact that the pure Newtonian contribution of the solvent is 
very important in the dynamic behavior of the CPyCl/NaSal solution. 

As for the compressibility, from Figure 2 the highest velocity in the fluid v m is approximately equal 
to 0.3m/s. Then, since the speed of sound in liquids is of order v s ~ 10 3 m/s, we have for the Mach 
number M ~ 10~ 4 . At such low Mach numbers, compressibility effects must be weak (PQ,[IH]), so the 
compressibility does not seem to be an important factor in the decay of the peaks. Let us present in the 
next section a detailed derivation that confirms this assertion. 



2 Actually, this is the same equation obtained for the Fourier transform of the velocity 



4 Influence of Compressibility 
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As mentioned above, at these low Mach numbers compressible effects must be weak. Nonetheless, in 
this section we'll be aimed at modeling the problem of sound waves propagation in viscoelastic fluids 
flowing inside pipes. The influence of viscoelastic properties in sound wave propagation phenomena in 
infinite fluids and in wall bounded flows (inside straight pipes) has already been addressed in |18j and 
[2], respectively. 

Let us first revisit the stress modeling problem in compressible fluids. It is well known (see for instance, 
PO] and [H]) that the motion in the neighborhood of any point consists of the superposition of: 

1- a uniform translation, 

2- a pure strain motion characterized by the rate of strain tensor D = Vi> -I- (Vw) , which can be 
decomposed itself into an isotropic expansion (contraction) in which the rate of extension of all line 
elements is §V-u and a straining motion without changes of volume characterized by — (|V • v) Sij, 
and 



3- a rigid body rotation with local angular velocity given by 



V x u 
2 



We shall restrict ourselves, as in previous sections, to flows in which the Weissenberg number is suffi- 
ciently small so we shall not consider shear-induced anisotropy in elastic stresses. Besides, we can assume 
a linear relationship between normal stresses and velocity derivatives by considering small oscillations 
of fluid particles in the wave motion. Moreover, as in the case of incompressible viscoelastic solutions, 
we shall consider that the stresses generated by the straining motion without changes in volume and by 
isotropic expansion relax on time scales Ai and A2, respectively. Then, the constitutive (linear) equations 
for the extra viscoelastic stresses may be written as follows 



and 



of + = - |m V • v I, (37) 

f)rr e 

al + X 2 -^-=r, 2 V-vI, (38) 
where I is the identity tensor of rank 2. The total viscoelastic stress a is obtained by summing <j\ and 

At this point, it is useful to note that for times larger than Ai the trace of the total viscoelastic stress 
satisfies the differential equation 

dTra 

Tra + X 2 ^^ = 3 m V-vI, (39) 

which coincides with the postulated equation for the trace of the stress tensor in Ref. [T5] when the 
displacements of fluid particles in the wave are small compared with the wave length or equivalently, 
when the velocity of fluid particles in the wave is small compared with the velocity of sound. 
The formal solution for the total viscoelastic stress can be written as 

/>oo 

a= / ds[G(s)I>(t- s) + {K(s)-l/3G(s)}ITrD(t- s)}, (40) 
Jo 

in which 

G(t) = (WAiK t/Al 



and 



K(t) = (r l2 /\ 2 )e- t ^ 



denote the relaxation functions of shear and normal stresses, respectively. 

Substituting this formal solution into the momentum balance equation leads to 

(v-\7)v+-^j=-\7p + T] s \7 2 v + J ds[G(s)\7 2 v + A(s)V(V • v)}, (41) 



where A(s) = K(s) + |G(s). 
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Since we are considering that the velocity of fluid particles is small compared with the velocity of 
sound we can neglect the term (v ■ V)v in the above equation. For the same reason, the relative changes 
in the fluid pressure and density should be small. Hence, we can write the latter magnitudes as 



P = Po+P / p = p + p', 

where p and p are the values of the pressure and density, respectively, in the unperturbed state, while 
the primed quantities refer to the oscillating part. 

The equation of continuity in the linear approximation can be written as 



p V ■v = 0. 



(42) 



We may also note, since we are considering small deviations from the state of equilibrium, that the 
changes in the entropy s are of second order of smallness |19j . In this approximation, to the first order 
of accuracy, we have that changes in the fluid pressure are only caused by variations in the fluid density. 
Then we can write p ~ p(p). Hence the small variations of p' are related to changes in p' in the following 
way 

'd£ 

dp 

where c is the speed of sound. 



P 



p'=c 2 P ', 



(43) 



After substituting for p' according to the above equation in 42 we obtain 



^ +c 2 poV -v = 0. 



(44) 



With the equations 44 and 41 for the four unknowns p' and v we are able to completely describe the 
sound wave. Moreover, in order to express the two unknowns in term of one of them it is convenient 



to take divergence in both sides of equation 41 and take the time derivative of equation [44] Then, it 
can readily be obtained the following equation describing the small oscillations of pressure (it was also 
derived in 18 ) 

o2 o poo Q 

K x V 2 p + Vs -V 2 p+ I dsM{s)-V 2 p(t-s), 



Po 



at 2 



at 



at 



(45) 



in which K x = p D c 2 is the compressibility modulus and M(s) = K{s) + |G(s). Here and henceforward 
we shall omit the prime symbol to the small variations of pressure. 

With the same assumptions as in section 3.1 concerning the flow geometry we shall consider that the 
velocity field in the sound wave can be described in terms of its component along the z axis, i. e. 



v = v z (r,z,t)e z 



(46) 



As a consequence of the latter, the pressure will vary only along the pipe axis and will be constant 
along a the cross section. 

p = p(z). (47) 

A noteworthy observation is that the velocity must depend on z, owing to the compressibility, which 
represent the main difference with the flow of an incompressible liquid. 

Let us seek a solution of equations [44] and [4l] in the form of a plane wave, 

v z (r,z,t) = Re[<j)(r)e likz - Wt) } p = Re[P e l{kz ~ wt) ]. (48) 
Substitution of the above expression for the pressure in equation [44] leads to the dispersion relation 



k 2 = 



2 r 



1 



l 



iwrjp frj s ^ \r) x 

poC x \r]o 3r] l — iw\ rj l — iw\ 



-i -l 



(49) 
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Similarly, after substitution of the expression for the velocity and the pressure itself in [41] we obtain 
an ordinary differential equation for <p(r) 



t^o 1 — icoXr]' iuj\rj' J 770 1 — iu>Xr) 




.kP 1-iuX 

-* — 1 — ^r^' ( 50 ) 



in which V? = ( r~ | , r? = 24 and ri' = 

r r d r y Qr J 1 1 '?« 12 w 

Solving the equation for (/>(r) we obtain for the velocity 



v z (r,z,t) = -- — 1- . (51) 



In the above expression it has been defined 



?7o 1 — iu>Xr}' iuj\r) 



Noting that 4 s = ikP e l ( kz wt \ we can conveniently rewrite the above equation as 

, .v 1 f, MM \ dp 

v z (r,z,t) = -- 2 —II— — , 52) 

™Po- HSx I ) + %) V MPcR)J dz 



which resembles the corresponding velocity for the incompressible case (Eq. 18). 
Hence, defining the local dynamic permeability as in section 3.2 by 

-r) v z (0,z,t) 



dp 

dz 



(53) 



and referring the frequency to J^, the inverse of the characteristic time of vorticity diffusion due to 
viscosity, it is obtained 

«c(«)=— y (54) 

1a - e c \ MPcjJ 

in which 

(■^A* \ia-n') + rh 

(55) 







2 


1(1 -»?') + % 




(1 - iaA) 






([f(W)+%] T^A + •/)) 



and 

^ = (56) 

In the above expressions it has been used Re x = . The dimensionless numbers A and a are defined 
in the same way as before. 

In experiments with the CPyCl/NaSal solution, the range of frequencies covered kept the number 
a in the range < a < 1. Estimating that the speed of sound in the solution is c x ~ 10 3 m/s, leads to 

Re x ~ 10 3 and thus (^rj ~ 10~ 6 . Let us also assume that ~ 1 (see [TH]) and that relaxation times 
of shear and normal stresses are of the same order of magnitude. Evaluating the remainder parameters as 
in the incompressible case we obtain that Abs(e c ) <C 1 leading to f3 2 = (3 2 . Thus, we reach the conclusion 
that for weakly compressible flows the LDR matches almost perfectly the corresponding incompressible 
flow. We can see very slight differences in the values of frequencies at which peaks are observed as well 
as in the amplitude of the LDR provided only that ~ 1 [2]. 



5 Conclusions 
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In this paper we have discussed some non-Newtonian effects of an Oldroyd-B fluid under oscillating 
pressure gradient. A relation was derived between flux velocity and pressure gradient in frequency 
domain that resembles Darcy's Law. An analytical expression was given for the dynamic permeability 
in terms of dimensionless parameters, which are: the number a that represents the ratio between the 
inertia and viscosity forces, the inverse of Deborah's number, which is the ratio between characteristic 
times of viscous effects and elastic ones, and the number wA — , which determines the differences between 
the Oldroyd-B and Maxwell models. Moreover, our calculations have predicted a transition from a Non- 
Newtonian (viscoelastic) to purely Newtonian behavior in the dynamics of viscoelastic solutions when 
the condition wA^ S> 1 is satisfied. This transition is therefore independent of the flow geometry, 
i.e. it has a universal character. We have found good agreement between the linear Oldroyd-B model 
predictions for the local dynamic response, as well as for the root mean square velocity at the pipe 
center, and the corresponding experimental data, shown in Figure 1 and 2, respectively. This suggests 
that it is the dissipative factor, present in Oldroyd-B fluids, which might leads to the decrease of the 
local dynamic response observed in experiments [3]. However, these good agreements are obtained when 
the solvent viscosity used in the fitting is considered as an increased coefficient owing to the possible 
fact that due to high concentrations of the viscoelastic solution, disturbances of the solvent velocity field 
generated by the motion of the polymer molecule (coil) affects the drag force on the neighboring coils. 
Consequently, the solvent viscosity is used here as a free parameter actually. The very small Mach numbers 
at which experiments were made indicate that compressibility effects can be neglected. Our derivations 
confirmed that, given the experimental conditions, weakly compressible effects won't contribute to any 
significantly change in the local dynamic response of the fluid compared against the incompressible result. 
However, discrepancies are evident at low frequencies in both figures. These disagreements in this range 
of frequencies could be produced by nonlinearities as the perturbation procedure developed suggested. 
According to our results, anisotropy in viscoelastic stresses should be stronger at low frequencies. The 
latter would drive the viscoelastic pipe flow weakly turbulent (through a nonlinear subcritical transition) 
if the fluid were perturbed significantly and consequently, in such state our hypotheses will break down. 
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